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Abstract 

The stability of the equation of state predicted by Skyrme-type interactions is 
examined. We consider simultaneously symmetric nuclear matter and pure neutron 
matter. The stability is defined by the inequalities that the Landau parameters 
must satisfy simultaneously. A systematic study is carried out to define interaction 
parameter domains where the inequalities are fulfilled. It is found that there is 
always a critical density p cr beyond which the system becomes unstable. The results 
indicate in which parameter regions one can find effective forces to describe correctly 
finite nuclei and give at the same time a stable equation of state up to densities of 
3 — 4 times the saturation density of symmetric nuclear matter. 
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1 Introduction 



Nucleon-nucleon effective interactions with an explicit density dependence have been 
largely employed since the seventies for studies of nuclear properties. Once the 
form of the interaction is chosen, either a zero-range force of Skyrme-type [Q] or a 
finite-range force of Gogny-type 0|, the parameters are determined by fitting some 
selected properties of doubly-magic nuclei and symmetric nuclear matter at satu- 
ration. Within a Hartree-Fock (HF) scheme these interactions are able to describe 
quantitatively the properties of stable as well as unstable exotic nuclei. The interac- 
tions are thus well controlled around the saturation density po of symmetric nuclear 
matter, for moderate isospin asymmetries and zero temperature. This type of ef- 
fective interactions has also been employed to study nuclear matter in conditions 
of astrophysical interest, for instance neutron matter with a finite proton fraction 
at densities up to several times po an d finite temperatures. This system at such 
conditions is relevant for studying proto- neutron stars. The effective interactions 
are thus extrapolated to conditions of density and isospin asymmetry which are not 
experimentally accessible, and one should then ask for the limits of validity of such 
an extrapolation. 

Of course, the description of nuclear matter consisting of nucleons is no longer 
valid beyond some value of the density, as other degrees of freedom appear. For 
instance, strange matter is expected beyond 3 — 4po ||. For the present discussion, 
we shall arbitrarily accept that for densities up to 4po, only nucleons are needed to 
describe nuclear matter. Even in this range of densities the effective interactions are 
not well determined. For instance, different effective interactions can give similar 
equations of state for symmetric nuclear matter and very different results for pure 
neutron matter. Indeed, the extreme asymmetry of isospin is not part of the usual 
input to determine the interaction parameters. In some cases, the determination 
of the parameters includes a fit to the equation of state of neutron matter calcu- 
lated with microscopic methods and realistic interactions. The Skyrme interaction 
RATP takes into account the variational microscopic results of Ref. ||. More 
recently, the set of SLy interactions || f?j have included the variational results of 
Ref. || among the conditions required to fix the parameters. It is worth noting 
that these microscopic variational calculations are not exempt from uncertainties 
since realistic interactions are not completely known, especially their three-body 
part. Indeed, the same theoretical model Q employing different realistic interac- 
tions leads to different equations of state, both in symmetric nuclear matter and 
in pure neutron matter, and sizeable differences appear for densities beyond the 
saturation value pq. Although this type of approach gives a precious guide for the 
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determination of the effective interaction parameters, it cannot replace the empirical 
data which are the natural input for the phenomenological interactions. 

The extrapolation of effective interactions can result in an unphysical behavior 
of nuclear matter [Sj. For instance, most Skyrme parametrizations predict that 
the isospin asymmetry energy e/ becomes negative when the density is increased. 
Consequently, the symmetric system would be unstable beyond the saturation den- 
sity, preferring a largely asymmetric system made by an excess of either protons or 
neutrons. Another type of instability refers to the magnetic properties of neutron 
matter. The possibility of a ferromagnetic transition at high densities has been 
studied long ago, employing different theoretical ingredients, and the results were 
contradictory. Currently used Skyrme interactions predict that neutron matter be- 
comes spin polarized at densities between ~ 1.1 — 3.5po H 0]- This transition 
would have important consequences for the evolution of a proto-neutron star since 
the mean free path of neutrinos would be zero []|, 11]. However, Gogny-type ef- 
fective interactions either exclude such a ferromagnetic transition or predict it at 
very high densities ||] . Relativistic mean- field calculations [12, [l3|, O] predict that 
such a transition could appear at densities beyond ~ 4po- Finally, recent Monte 
Carlo simulations [15] as well as Brucckncr-Hartree-Fock calculations [JH], 17] using 
modern two- and three-body realistic interactions do exclude such an instability. 

The purpose of this work is to analyze whether the presence of instabilities is 
inherent or not to the Skyrme- type interactions. In Section 2 we explain the restric- 
tions on the Skyrme parameters imposed by the requirement that no instabilities 
should appear at densities in the range (1 — 4)po- Furthermore, the sound velocity 
v s /c in matter must remain smaller than unity, and this adds a new constraint. 
Section 3 contains the discussion of the results. Conclusions are drawn in Section 
4. Useful expressions are given in Appendix A and B. 



2 Constraints from the Landau parameters 

The condition that the spin-unpolarized neutron matter should be the most stable 



phase at any density was used in Ref. [18] to constrain the Skyrme parameters. As 



the potential energy of the fully spin-polarized phase only depends on the combi- 
nation £2(1 + X2), ferromagnetic collapse is necessarily avoided if this combination 
is positive. Taking into account current Skyrme parametrizations, the authors of 



Ref. [18] concluded that —5/4 < X2 < —1. However, this condition is not sufficient: 
the family of SLy parametrizations || [jj imposes this constraint but at densities of 
about 0.37 fm -3 , a ferromagnetism instability is predicted in neutron matter |]]. 
Here, we make a systematic use of the stability criteria related to the adimen- 
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sional Landau parameters Fi, Gi, F t and in symmetric nuclear matter, and iy 

(n) 

and G\ in neutron matter. These criteria simply establish that any parameter of 
multipolar ity I should be greater than —{21 + 1). Current effective interactions sat- 
isfy these conditions at saturation, and we shall investigate if they can be maintained 
for densities up to 4po- 

The general form of the Skyrme interaction contains ten parameters, eight of 
them denoted as ti, Xj (with i = 0,1,2,3), the other two being the power a of 
the density dependence and the strength Wq of the zero-range spin-orbit term. 
These parameters are usually determined by fitting some selected properties of finite 
nuclei and symmetric nuclear matter. In this study we follow a different procedure, 
using some accepted nuclear matter values to fix as many interaction parameters as 
possible and using the stability criteria to put bounds on the remaining parameters. 

Our starting point is given by the empirical values of the following quantities in 
symmetric nuclear matter at the saturation point: density po, energy per particle eo, 
effective mass m^, compression modulus Kq, surface energy eg, and isospin asym- 
metry energy ej. Imposing that a Skyrme- type interaction reproduces the values of 
po, eo, m,Q, and Kq completely determines the parameters to, t^, a and the combi- 
nation To = [3*i + (5 + 4x2)t2]/8. The explicit expressions are given in Appendix A. 
With a Skyrme interaction, the surface energy es depends on the previously fixed 
parameters, on the combination T$ = [9t\ — (5 + 4x2)*2]/8 and on the strength Wo 
of the spin-orbit term [jOJ. The value of Wo is usually determined by fitting the 
spin-orbit splitting of some selected levels in finite nuclei, and for the present dis- 
cussion we shall assume the often used value Wo = 120 MeV-fm 5 . Furthermore, the 
symmetry energy ej can be expressed in terms of the above parameters and of the 
four combinations ijXj (i = — 3). In summary, from the ten Skyrme parameters, 
one is kept fixed (Wo), six parameters or combinations are determined by the six 
empirical inputs, and three combinations are free. We found convenient to choose 
the free combinations x = t\X\, y = £2^2 , z = £3X3. 

From here on, our task is to analyze the possible domains of the (x, y, z)-space 
for different densities. For a given value of the density p we would like to know 
what is the (x, y, z) domain where no instability can occur. Stability implies that 
any adimensional Landau parameter of multipolarity I must be larger that — (2/ + 1). 
Skyrme forces only contain monopolar and dipolar contributions to the particle-hole 
interaction so that all Landau parameters are zero for I > 1. Thus, we have twelve 
inequalities, eight coming from symmetric nuclear matter, and four from neutron 
matter. Explicit expressions of the twelve Landau parameters are given in Appendix 
B. 

Two of the twelve inequalities, however, play a special role. The condition 
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F\ > —3 is trivially satisfied for any density as long as < m*(po) < m, which is 
the case. The Landau parameter To only depends on t%, t%, a and the combination 
To, which are determined by the initial inputs. Exploring the inequality Fo > — 1 
as a function of the density, it is found that it is violated at a density smaller than 
po. It corresponds to a density where the compression modulus becomes negative, 
which is the spinodal point, i.e, the occurrence of a liquid-gas transition. Note 
that within a Skyrme interaction framework, the spinodal density is completely 
determined by the saturation properties eo, po, Kq and mg. After excluding the 
inequalities related to F\ and To, we are left with ten inequalities to be satisfied by 
the yet free combinations x, y, z. 



2.1 Symmetric nuclear matter 

The parameters G' and G' x give two constraints for y: 

10d(p) 2 m 

10Co(p) 2 

y > + -(T -2Ts , (2 

p 3 

where the Cj and «j are defined in Appendix A. The parameters F{, G\ give two 
constraints on combinations of x and y: 

^-7^ > + — (T -2T 5 , (3 

5 p 15 

3 4C (p) 4 

z + < — T -2T 5 ). (4 

5 p 15 

The remaining Landau parameters Go and Tq constrain combinations of all three 
x, y and z unknowns: 

±-{p° - pl)z + (p 2 / 3 - )x + \( P ^ + P T)y > (5) 

^_1^2/3 , 2/3, 



— (c.ip) + d(po) + ^)~ y& + p ° )(T ° - 2Ts) ' 
- po> - (p 2/3 - po /3 )- + \(p 2 " - pT)v > (6) 
-i- (cm - Cl { PQ ) - *L) - Up^ _ p f )(To _ 2Ts) . 

3ai \ po J 15 

However, it is possible to take the linear combination Tg + Go to get another con- 
straint independent of z. Note that Eq. (||) reduces to the trivial condition ej > 
when p = pq. It is due to the fact that the symmetry energy e/ has been used to 
write xq in terms of the remaining parameters. 
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2.2 Neutron matter 

The parameters F^ n \ and give two constraints on combinations of x and y: 

X ~h < ^ + ^o-2T s ), (7) 

' < ?4> _A(T _ 2rs) . (8 , 
5 p 15 

The remaining Landau parameters i 7 ^, and Gq™^ constrain combinations of all 
three x, y, z parameters: 

(, + l)(, + 2 K-M 3 
6(4a 2/ o 2 / 3 - 3aiPo /3 ) 5 



4a 2/ o 2 / 3 - 3aip 



2/3 



Gi(po) - G 2 (p) + ^) + -i(T - 2T S ) , 



±(fT - Po > + ^(2a 2 p 2 / 3 - 3ai^ /3 )x + ^(2a 2/0 2 / 3 + 3a iP 2 /3 )y > (10) 
4 'dipo) + G 3 (p) + ^) - ^(7a 2 p 2 / 3 + 3a lP 2/3 )(T - 2T S ) . 



3 V v ™ vr/ Po/ 45 

Note that the linear combination cki-Fq + ^0™^ gi ves another constraint independent 
of z. 

2.3 Sound velocity constraint 

In addition to the constraints from stability requirements it is important to check 
that for each density the sound velocity v s remains smaller than the speed of light, 
i.e., superluminosity does not occur. The sound velocity is directly related to the 
compression modulus K(p) which can itself be expressed in terms of the Landau 
parameters Fq and F\\ 

mv 2 = -K (11) 
9 v ; 

h 2 k 2 F 1 + F 
3m 1 + ±Fi ' 

In symmetric nuclear matter it turns out that v s depends only on the parameters 
to, £3, <y and the combination To, all of which being already determined by the six 
initial empirical inputs. Therefore, there will be no additional bounds on the allowed 



volume in (x, y, z)-space. However, Eq. (11) shows that v s increases with increasing 



density and for a given choice of initial empirical inputs there is always a value 
Psound beyond which v s /c is greater than 1. The value of p soun d depends essentially 
on the adopted value for Kq. For instance, p sound is 

around 3.5p if K = 350 MeV 
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and it may become 6po if Kq = 250 MeV. Thus, there is no Skyrme interaction 
which can be reasonably be used beyond p SO und because it would predict unphysical 
values of the sound velocity. 

The situation is somewhat different in pure neutron matter. Now, the Landau 
parameters in Eq. ([ll]) depend also on the (x, y, z) parameter combinations. The 
requirement that v s /c remains smaller than unity adds one more constraint to the 
determination of the allowed volume in the (x, y, z) space. We shall see in the 
next section that the sound velocity in neutron matter does not bring any effective 
restriction on the allowed volume in the parameter space, at least around the values 
of the six empirical inputs used in the present study. 

We must mention that a study of the sound velocity in asymmetric nuclear 
matter was made in Ref. pG|]. The aim was to find the conditions to be satisfied by 
the Skyrme parameters so that superluminosity would never appear at any density 
whatsoever. The conditions which were found are very restrictive. In the present 
work, our point of view is quite different since we don't think that the Skyrme 
effective approach should be valid at densities beyond 4po- 

3 Results and discussion 

The inputs used in the present study are: eo = —16.0 MeV, po = 0.16 fm~ 3 , 
Kq = 230 MeV, m* /m = 0.70, e s = 18.0 MeV, and ej = 32.0 MeV. 

For a given value of the density p, the various constraints define in the (x, y, z) 
parameter space an allowed volume that we call Q. The surface of f2(p) is made of 
planes because all constraints are linear in x, y, z. Any point outside f2 is forbidden 
because some of the inequalities would not be fulfilled. To present the results and 
facilitate the discussion we will explore the z-axis. The intersection of a volume 
O(p) with a z = constant horizontal plane gives a polygon which can be represented 
in the (x, y)-plane. By varying p in the range p$ — 4po one can follow the evolution of 
the polygons. The vanishing of the area of the polygon at some critical density p cr 
indicates that there is no Skyrme interaction (having the chosen value of z) which 
can fulfill the chosen set of constraints at densities beyond p cr . 

We have performed two types of calculations: a) using only the six constraints 
from the Landau parameters of symmetric nuclear matter and disregarding those 
of neutron matter; b) using all eleven constraints from symmetric nuclear matter 
and neutron matter plus the sound velocity constraint in neutron matter. In Fig. 1 
are displayed some typical results. We have chosen a positive and a negative value 
of z since there is no a priori limitation on z. The value z = 2 • 10 4 MeV fm 3 ( 1+CT ) 
correspond to Figs, la) and lb), and z = — 2 • 10 4 to Figs, la') and lb'). Results 
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involving only Landau parameters of symmetric matter are displayed in Figs, la) and 
la'), whereas Figs, lb) and lb') correspond to results using all eleven constraints. 
The contours are drawn every 0.5po starting from pq. 

Let us first examine the results without neutron matter constraints (Figs, la 
and la'). At p = po the contour is a polygon whose upper and lower horizontal 
sides are determined by the constraints on G' and Go, respectively, whereas the left 
and right sides correspond to F[ and G\. When p increases the lower side becomes 
tilted and the constraint on Fq gradually appears as a new side on the right hand 
side of the polygon. At p ~ 2po and above the Fq constraint is dominating over 
the G\ constraint. The surface of the polygon shrinks as p increases. The value of 
p cr is above 4.5po- This value is reached for x and y both negative if z is positive, 
whereas for z negative p cr corresponds to x > and y < 0. 

The situation changes somewhat when neutron matter constraints are added 
(Figs, lb and lb'). At p = po the lower side is given by the Gq constraint while 
the right side of the polygon is mainly determined by the Fq constraint which 
limits severely the allowed area. It must be noted that the sound velocity does not 
bring any limitation on the (x, y, z) parameters in the domain explored here. One 
can see that the critical density p cr has now a lower value as compared to the case 
of Figs, la-la'. 

There is no natural limitation to the domain of the z parameter and therefore, we 
must explore the dependence of the results on z. In Fig. 2 we present calculations 
where all eleven constraints are included, exploring the z-axis from negative to 
positive values. The interpretation of the different sides delimitating the contours is 
the same as in the case of Fig. 1. One can see that p cr is highest when z is between 
and 3 • 10 4 . For negative values of z the allowed area at p = 2.5po is already fairly 
small. Below z — 3 TO 4 the polygons of larger p are contained inside those of smaller 
p, i.e., interactions stable at density p are also stable at all densities between po and 
p. For z larger than 3 • 10 4 the allowed areas at larger p tend to move outside the 
area of p = po, which means that that such interactions might seem acceptable at 
large densities but they have the defect of having instabilities at normal density. 

So far we have discussed the results obtained with the values of empirical inputs 
as adopted in Section 2. It is interesting to see how the conclusions may depend 
on this choice. In Fig. 3 we study, for a fixed value of p = 2po and a chosen 
z = 10 4 MeV.fm 3 ( 1+<T ) , the evolution of the allowed polygons when one changes one 
of the empirical inputs: the effective mass ra^jm (Fig. 3a), the incompressibility K$ 
(Fig. 3b), the surface energy €5 (Fig. 3c). It can be seen that the results are sensitive 
to the value of mg/m. A smaller value of iri^/m makes the allowed domain larger 
and therefore, the critical density p cr is shifted to higher values. On the contrary, 
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the allowed domain is reduced when higher values of m^/m are considered, and 
instabilities appear for densities below 2pQ. Hence, to avoid low density instabilities, 
one should prefer low mg/m. The sensitivity to eg is also non negligible, whereas 
the dependence on Kq is moderate. 

In Fig. 4 we present the calculated values of the critical density p cr (in units of 
po). Fig. 4a shows p cr as a function of the parameter z for various choices of the 
empirical input mg/m around the value 0.70 . As pointed out above, one can see 
clearly that it is more favorable to choose z positive and less than 3 • 10 4 in order 
to have p cr not too small, and that larger values of mj/m tend to yield lower p cr . 
The three curves of Fig. 4b represent p cr as a function of the relative variations 
(in percentage) around the standard choice of mj/m, Kq and eg- The value of z is 
fixed at 10 4 MeV.fm 3 ( 1+,J ). These curves show that p cr depends little on Kq whereas 
moderate increases of mg/m (from 0.70 to 0.75), or of eg (from 18 MeV to 20 MeV) 
can lower p cr below 3pQ. 

Now, the question arises whether these bounds for Skyrme parameters are use- 
ful or not for calculations in pure nucleon matter and finite nuclei. In principle, we 
expect a qualitative positive answer for finite nuclei because the used inputs guar- 
antee that the first terms in a leptodermous expansion of the mass formula are well 
described. To be more quantitative, we have constructed several parametrizations 
using the following criteria. We first fix an arbitrary value of z, going from —10 to 
4- 10 4 , and then we take the values of x and y as the coordinates of the point Q(p cr ). 
In this way instabilities are pushed to the highest values of the density. The sets of 
parameters are dubbed Skzn and they are given in Table |]. Of course, the criteria 
of highest value of p cr is purely arbitrary but our intention here is just to explore 
the resulting parametrizations and not to give the best parameter set. In Fig. 5 are 
displayed the results for pure neutron matter. One can see that choosing negative 
values of z results in an increasing neutron effective mass for increasing densities. 
On the other side, values of z greater than 2 • 10 4 produce too low binding energies. 
We have then restricted our calculations in finite nuclei to parametrizations SkzO, 
Skzl and Skz2. In Table § are displayed the ground state binding energies and 
charge radii of some doubly-magic nuclei together with the results obtained with 
SIII and SLy4 interactions. The results are reasonably satisfactory. We have no- 
ticed that convergence problems occur with SkzO when looking for the HF solution 
of 208 Pb. To obtain better results one should relax the condition of highest value 
for the critical density, and of course include some finite nuclei results in the fine 
tuning of the parameters. 
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4 Conclusion 



In this work we have explored the parameter domains of Skyrme-type forces where 
the stability conditions related to Landau parameters inequalities can be satisfied 
up to densities of the order of 4po- Stability of both nuclear matter and pure neutron 
matter are considered. We have taken advantage of the possibility to characterize 
the domains by analytical expressions. Starting from a general Skyrme force con- 
taining 10 parameters, we have shown that 7 parameters or combinations thereof 
can be approximately fixed by physical quantities which can be considered experi- 
mentally known to some extent. Thus, the problem reduces to the study of allowed 
domains in a 3-dimensional space spanned by 3 parameter combinations. 

The results show that, for any Skyrme-type interaction there is a critical density 
p cr above which one cannot insure all stability conditions. This critical density does 
not exceed 3.5 — 4/?o for a reasonable choice of empirical inputs. It is nevertheless 
an interesting result to know that it is possible to find Skyrme-type interactions 
which can give stable nuclear matter and neutron matter up to such densities. The 
parameter domains are well identified and it would be worthwhile to look inside 
those domains for Skyrme interactions which can also describe accurately finite 
nuclei. Our exploration with parametrizations Skzn has shown that reasonable 
ground state binding energies and charge radii can be obtained. A more systematic 
study over a wide number of nuclei would be necessary to get the Skyrme-type 
interaction which is stable over the largest range of densities. This new interaction 
would be very useful for neutron star calculations, like for instance the neutrino 
mean free path or the URCA process. 

The present analysis is restricted to Skyrme interactions, because the simplicity 
of the resulting expressions allows for an algebraic study. The presence of instabil- 
ities in nuclear matter and neutron matter beyond some critical density could be 
attributed to the zero-range form of the interaction, and its /c 2 -dependence of the 
effective masses. Finite range effective interactions Q are an alternative choice for 
describing the equation of state of nuclear and neutron matter. We know that the 
actual parametrizations (Dl, D1S, DIP) do not predict instabilities || for densities 
below ~ 5po. We have not considered a systematic study of Gogny effective inter- 
actions because they contain more parameters, and they do not allow for a simple 
algebraic analysis. 
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A Nuclear matter relations 

Employing a Skyrme interaction, the energy density of semi-infinite symmetric nu- 
clear matter is written as 



H 



n 



i 



m 



2m^ + 8^ + 16^ + 8^' " W V - PW 



(12) 



where V so = 9W 2 /16. 

Denoting by po, €q, K and respectively the density, energy per particle, com- 
pressibility, and effective mass at saturation of the bulk symmetric nuclear matter, 
it turns out that the Skyrme parameters to, t$, a, and the combination To can be 
written as: 



2*1 _?n i 

ml m J po 



(13) 



-K + e + 



h 2 



2ti 2 



10m 15m 



4(0) 



a 



o , 



-e + 



3ti 2 



10m 5m, 



4(o) 



o , 



(14) 



16 1 



-e + 



3h 2 



h 2 



10m 5m, 



4(o) 



o , 



(15) 



to = 

3po 



3h< 



eo 



*4(0) - T^aPo 



l+cr 



Wm* r " ' 16 



(16) 



where po = 2/(3tt) 2 kp(0). Once the parameters to,t^,a and the combination To 
have been fixed, the surface energy only depends on the combination Ts and the 
spin-orbit strength Wo- The surface energy can be written as an integral over the 



density |19| 



87rr, 



po 

) / d P 

o 



Tr 



5 T oP + \ T sP 



36m 36 
^ + 8^ + 16^ 



•mr 2 

w v soP 



1/2 



(17) 



1/2 
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where ro = [3/(47rpo)] 1 ^ 3 is the unit radius, and a Thomas-Fermi approximation up 
to ?i 2 -order has been used to replace the kinetic energy density. 

Then, the following relation involving the symmetry energy e/ can be used to 
relate xq to the already known parameter combinations and to the free parameters 
(x,y,z): 

t x = + + ^(T - 2T 5 )) a lP 2 /3 - ± Zf % - 2C 1 ( /9o ) - ^ • (18) 

To simplify the notation used in Section 2 we have introduced the following 
functions of the density: 

Co(p) = — + T p, (19) 

771 

Ci(p) = -(^ + T op)^ + \to + ^ , (20) 



chip) = -\{^+^Ajr^-\^-k {(T + l){a + 2)hpcJ (21) 



~ / x 1 I ^ \ "2 1 1 

CM = -2^ + r ^J^73 + 2 t0+ 12^ ' (22) 

which contain only known parameters. We have also introduced the constants 

/ 4 \ 1/3 / 4 \ 1/3 

« = ■ <-(t) < 23 » 

B Landau parameters 

The Landau parameters in symmetric nuclear matter are: 

„ /3 1 , . , „\ 2m* kr? „ 

F ° = (4*o+ig(^+l)(^ + 2)*3P a )-^f-Fi , (24) 

/l 1 \ 2m* kw 

G = ^-t (2x -l) + -t 3 p^2x 3 -l)j^-f -Gi , (25) 

= ("^0(2x0 + 1) - -i 3 p <T (2x 3 + 1)J -^f , (26) 
1 . 1 . „\ 2m*k F 



Co = ( -4*0 " ^ ) -^"Gi. ( 27 ) 
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Fx 
G x 
F[ 

where the Tj's are the following parameter combinations: 

T = ^ [3*i + 5t 2 + 4y] , 

T x = \[2x + 2y-tx+t 2 ] , 
o 

T 2 = -[2x-2y + tx-t 2 ] , 
o 

t 3 = lh-t 2 ] , 

and the effective mass is: 

— = - + T oP . 

m* m 

The Landau parameters in neutron matter are: 
F (n) = (lt (l-x ) + j l (a + l)(<j + 2)t 3 p°(l-x i 
G { o n) = Q^o-l) + ^3P CT (x 3 -l))^-G 

F (n) = _ 3(To _ T2) < P) 

n 

Gt ] = -3(T 1 -r 3 )^p, 
and the effective mass is: 

— = — + {T - T 2 )p . 
m* m 



= -3T — g-p , 
n 

m 

-3Tx W p , 

* 

m 

= 3T 2w p, 
m 

= m^p, 
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z 

MeV fm 3 ( 1+CT ) 


X 

Mev fm 5 


y 

MeV fm 5 


*2 

MeV fm 5 


X 


X\ 


%2 


x 3 


Skz-1 


-10 4 


570.42 


266.2 


-299.14 


-0.2665 


1.2968 


-0.8899 


-0.7282 


SkzO 





458.0 


215.0 


-258.18 


0.1986 


1.0413 


-0.8328 


0.0 


Skzl 


10 4 


217.0 


221.0 


-262.98 


0.6052 


0.4933 


-0.8404 


0.7282 


Skz2 


2 ■ 10 4 


24.0 


227.0 


-267.78 


1.0290 


0.0546 


-0.8478 


1.4564 


Skz3 


3- 10 4 


-265.5 


233.1 


-272.66 


1.4174 


-0.6082 


-0.8549 


2.1846 


Skz4 


4- 10 4 


-512.0 


239.15 


-277.50 


1.8226 


-1.1640 


-0.8618 


2.9127 



Table 1: Skzn parametrizations. All of them have the following common values: a=0.1694, 
to = -2471.10 MeV fm 3 , h = 439.85 MeV fm 5 , t 3 = 13732.8 MeV fm 3 ^ 1+a \ as they come 
from the initial inputs. 



Nucleus 


Observable 


SkzO 


Skzl 


Skz2 


Exp. 


SIII 


SLy4 


16Q 


B/A(MeV) 


8.39 


8.38 


8.38 


7.98 


7.95 


7.97 




r c (fm) 


2.76 


2.76 


2.76 


2.73 


2.76 


2.80 


40 Ca 


B/A(MeV) 


8.87 


8.85 


8.85 


8.55 


8.48 


8.55 




r c (fm) 


3.47 


3.47 


3.47 


3.49 


3.50 


3.51 


208 pk 


B/A(MeV) 


N.C. 


7.97 


7.91 


7.87 


7.79 


7.80 




r c (fm) 


N.C. 


5.48 


5.49 


5.50 


5.58 


5.52 



Table 2: HF binding energies and charge radii for several Skzn parametrizations. N.C. 
means that the HF calculation has not converged. 
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Figure 1: Comparison of the sections of the volume fi(p) by a horizontal plane z = 
constant, in units of MeV fm 3(1+CT ) . The horizontal and vertical axis are for the parameters 
x and y, respectively, in units of MeV fm 5 . Cases a) and a') correspond to the constraints 
from the Landau parameters of symmetric nuclear matter. Cases b) and b') include also 
the constraints from neutron matter. Two values of z have been used, namely z = 2 ■ 10 4 
for ) and b), and z = —2 • 10 4 for ') and b'). The different closed contours 

correspond to different values of p. The largest area is for p = p , the next one corresponds 
to an increase by a step of 0.5pcb an d so on. For the sake of clarity, they are alternatively 
represented by solid and dashed lines. 
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z = - 1.5 10 4 z = z=1.5 10 4 




Figure 2: Comparison of the sections of the volume fl(p) by a horizontal plane, for choices 
of z different from Fig. 1. The various closed contours correspond to different values of p. 
The largest area is for p = p , the next one corresponds to an increase by a step of 0.5po, 
and so on. The calculations are done using the full constraints. The values of z are in 
units of MeV.fm 3 ( 1+,J ). 
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Figure 3: Sensitivity to some empirical inputs. Case a: results for mj/m = 0.6, 0.7 

and 0.8 (respectively dotted, solid, dashed), keeping the remaining empirical inputs to 

their fixed values. Case b: same as a), with K = 210, 230 and 250 MeV (respectively 

dotted, solid, dashed). Case c: same as a), with the surface energy €$= 16, 18 and 20 

MeV (respectively dotted, solid, dashed). All cases have been calculated for z = 10 4 MeV 
fm 3(i+<x) and p = 2po _ 
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Figure 4: Case db'. Per ^ 

function of z (in units of MeV.fm 3 (l + a)) for different values 
of mQ/m=0.6,0.7,0.8 (respectively dotted, solid, dashed). Case b 

• Per ^ 

function of 

relative variations of some empirical inputs; dotted line: variation of m^/m, solid line: 
variation of K , dashed line: variation of eg. 
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Figure 5: Predictions of several Skzn parametrizations for the binding energy and the 
effective mass of pure neutron matter and for asymmetry energy as function of the density 
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